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We present our recent semi-analytical evaluation of the pion-pole contribution to 
hadronic light-by-light scattering, using a description of the pion-photon-photon 
form factor based on large- No and short-distance properties of QCD. Inclusion of 
all light pseudoscalar states leads to a bbyLiPS = +8.3 (1.2) x 1CT 10 . We also sketch 
an effective field theory approach to hadronic light-by-light scattering. It yields 
the leading logarithmic terms that are enhanced by a factor No and shows that 
the sign of a„ ? 1 cannot be inferred from a constituent quark-loop. In view of 
several problematic issues that still remain to be clarified, our estimate for the full 
hadronic light-by-light scattering contribution is a^ byL,had = +8 (4) X 10 — 10 . 



1 Introduction 

In early 2001, the Muon g — 2 Collaboration at the Brookhaven National 
Laboratory announced u a deviation of 2.6 a between their measured value 
for a M from the prediction in the Standard Model (SM). 

At the same time, we were working on a description of the pion-photon- 
photon transition form factor JF 7r o 7 » 7 . based on large-A^ QCD, which in- 
corporates short-distance constraints from the operator product expansion 
(OPE) 13. This form factor enters in the numerically dominant pion-pole con- 
tribution to the. hadronic light- by-light scattering-correction to g — 2, denoted 
by aJ; byL;had . una Since the form factors used in Era did not fulfill all short- 
distance constraints, we found it worthwhile to reevaluate this pion-pole con- 
tribution a^ byLi7r °. Furthermore, .we wanted to push the analytical calculation 
further than was done in Refs. Era, which relied purely on numerical methods. 

Our semi-analytical calculation El produced a result with roughly the same 
absolute size-bjji. with an opposite sign compared to the evaluations by three 
other groups LI'dQ. In an accompanying paper we therefore provided a deriva- 



*Talk presented at the Workshop on the Phenomenology of Large-Ap QCD, Tempe, 
Arizona, USA, 9—11 January 2002. To be published in the Proceedings. 



March 26, 2002 



tion of the coefficient of the leading log-square term, based on an effective field 
theory (EFT) approach and the renormalization group (RG), which lead again 
to a positive coefficient. The correctness^ the positive result for a^ byL;7r ° was 
later confirmed by the authors of Refs. uu and by independent calculations, 
see Ref. U. The new result for a^ yh ' haji reduces the discrepancy between the 
experimental value and the SM predicticoi to about 1.6 <j. 

Below I will present our evaluation u of a^ byL;7r ° (Sec. ||) and sketch an 
EFT approach to a ^y^ aA (Sec. ^). Section [| contains some considerations 
about leading and next-to-leading logarithms and an argument why the sign 
of a LbyL;had canno t b e inferred from a constituent quark-loop. I will conclude 
with an estimate for a ^y^^ d with a conservative error bound which takes into 
account the still unsolved problems. 



2 Pion-pole contribution 

The quantity a Lb y L ; had mvo i ves tfi e fourth-rank vacuum polarization tensor 
n At! yAp('Zi, 92, 93), built from the light-quark electromagnetic current = 
(2u7 M u — d-y^d — sj f _ t s)/3. In this section we will concentrate on the con- 
tribution from the neutral pion intermediate state, see Fig. |l|, which can be 




Figure 1. The pion-pole contributions to light-by-light scattering. The shaded blobs repre- 
sent the form factor jF 7r o 7 * 7 * {qf, defined in Eq. (^). 
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The pion-photon-photon transition form factor J r 7I .o 7 * 7 * is defined by 



(1) 
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and T\ and T2 in Eq. (|l]) are polynomials of up to sixth order in the momenta 
p, qi , and q 2 , with p 2 — m 2 , see Ref. □ for details. 

Since no data on the doubly off-shell form factor J r 7r o 7 » 7 . (q 2 , q%) is avail- 
able, one has to resort to theoretical descriptions of the form factor. In order 
to proceed with the analytical evaluation of the two-loop integrals in Eq. (|l|) , 
which is very difficult without specifying J r w o 7 » 7 « {q 2 , q 2 ), we considered a cer- 
tain class of form factors which includes the ones based on large- Nc QCD 
that we studied in Ref. □. For comparison, we have also used a vector meson 
dominance (VMD) form factor (sdiich leads to a result close to the more so- 
phisticated models used in Refs. ma), and a constant form factor, derived from 
the Wess-Zumino-Witten QVZW) term, that describes the chiral anomaly in 
chiral perturbation theory In large- Nc QCD, the pion-photon-photon form 
factor has the form 

-r- 1 2 2\| _ c »j(ffl; ffl) /q\ 

where the sum runs over an infinite set of narrow vector resonances. The 
polynomials Cij(q 2 ,q 2 ) are arbitrary, however, there are constraints at long 
and short distances. The normalization is given by the WZW term, 
J r 7r o 7 . 7 ,(0,0) = -Nc/QSi^Fn), whereas the OPE tells us that 



ft >>~rA*f. <*-*?) - I £ + £ *j? + o 

In the following, we consider the form factors that are obtained by truncation 
of the infinite sum (||) to one (lowest meson dominance, LMD), respectively, 
two (LMD+V), vector resonances per channel: 

l?1 ' q%) ~ 3 (q 2 - )(q 2 -M 2 )' [b) 

Ft MD + v^ _2x _ *V gi 2 gl(g?+g 2 2 ) + fei(g?+g 2 2 ) 2 + h 2 qhl + Mgi+g|) + h 7 

^o 7 , 7 ,l9l,9 2 J 3 (g2 _ M 2 i)(g2 _ M 2 2)((? 2_ M 2 i)(9 2_ M 2 2) > 

(6) 

with c y = N c M v /(4tt 2 F 2 ) ,h 7 = — N c M Vi My 2 /(4tt 2 F 2 ). Not all parame- 
ters in the LMD+V form factor are fixed by the normalization and the leading 
term in the OPE. We have therefore detexmined the coefficients hi, /i 2 , and h§ 
phenomenologically. According to Refs. E30 the form factor T^o^*^* (— Q 2 , 0) 
with one photon on shell behaves like l/Q 2 for large spacelike momenta, 
Q 2 = —q 2 . Whereas the LMD form factor does not have such a behavior, 
it can be reproduced with the LMD+V ansatz, provided that hi = 0. A 
fit of the LMD+V form factor to the CLEO data M, with hi = 0, yields 
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/15 = 6.93 ± 0.26 GeV 4 , see Ref. i. Analyzing the experimental data El for the 
decay tt° — ► e + e~ in the context of Ref. Oj, leads to \h 2 \ ^ 20 GeV 2 . 

The WZW form factor is given by . («i > ?I) 

the usual VMD form factor reads 



-iV c /(127T 2 i^) and 



fVMD / 2 2\ 



N c 



Ml 



Ml 



\W V {q\ - Ml) (qj 



Ml) 



(7) 



Note that this form factor does not correctly reproduce the OPE in Eq. 
All form factors above can be written in the following way 



3 



/(«?) - E 



1 



4 



9M Vi {ql) 



(8) 



where the explicit expressions for f(qf) and gu v . (q±) for the different form 
factors can be found in This crucial observation allows one to cancel all 
dependences on qi-q 2 in the numerators in J r 7r o 7 . 7 » (q 2 , ((Z1+92) 2 ) in Eq. (Q) and 
to perform all angular integrations in thc_two-loop integrals analytically using 
the method of Gegenbauer polynomials E-i The pion-exchange contribution to 
a p can then be written as a two-dimensional integral representation as follows, 
where the integration runs over the moduli of the Euclidean momenta: 



LbyL; 77 



(9) 
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(10) 



(I 



LbyL;,r°(2) 
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M=M 7r ,M Vi 



The dependence on the form factors is contained in the functions f^ 1 \g i M v ,, 
and g { M whereas the universal [for the class of form factors with the decompo- 
sition (^)] weight functions Wf 1 ,w gi ,Wg 2 are rational functions, square roots, 
and logarithms a. We have plotted them in Fig. ||. The functions and w gi 
are positive and concentrated around momenta of the order of 0.5 GeV. Note, 
however, the tail in in the Q\ direction for Q 2 ~ 0.2 GeV, which pro- 
duces for the constant WZW form factor a divergence of the form ~ C In 2 A 
for some UV-cutoff A. From the analytical expression for Wf 1 we obtain El 
C = 3[Vcm/(127ri 7 ' 7r )] 2 . On the other hand, the function uu g2 has positive and 
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Table 1. Results for the terms aJ; byL; " ^\ a L 1 byL, ' r ^ and a J; byL;,r0 according to Eq. |j) 
for the different form factors. In the WZW model we used a cutoff of 1 GeV in the first 
contribution, whereas the second term is ultraviolet finite. In the LMD+V ansatz we have 
set hi = GeV 2 and h 5 = 6.93 GeV 4 . 



LbyL;x U (l) 



Form factor 



a 



a 



LbyL:: 



a 



LbyL;- 



x 10 



10 



WZW 
VMD 
LMD 

LMD+V {h 2 = -10 GeV 2 ) 
LMD+V (h 2 = GeV 2 ) 
LMD+V {h 2 = 10 GeV 2 ) 



0.095 
0.044 
0.057 
0.049 
0.045 
0.041 



0.0020 
0.0013 
0.0014 
0.0013 
0.0013 
0.0013 



12.2 
5.6 
7.3 
6.3 
5.8 
5.3 



negative contributions in the low-energy region, which will lead to a strong 
cancellation in the corresponding integrals. 

In Table [l] we present the numerical results for the different form factors. 
One observes that all form factors (apart from the unrealistic constant WZW 
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form factor, which only serves illustrative purposes) lead to similar results, 
given mainly by a^ byL,,r ^ . It is important to correctly reproduce the slope 
of the form factor at the origin and the data at intermediate energies. On 
the other hand, the asymptotic behavior at large Qi seems not very relevant. 
All these observations can easily be understood from the plots of the weight 
functions. For a fixed value of /12 in the LMD+V form factor, our results are 
rather stable under the variation of the other parameters. If we vary Mv t , 
M V2 , and h 5 by ±20 MeV, ±25 MeV, and ±0.5 GeV 4 , respectively, the result 
for a^ byL; "° changes by ±0.2 x 10~ 10 . In contrast, if all other parameters are 
kept fixed, our result depends for \Ji2\ < 20 GeV 2 almost linearly on hi- In 
this range a^ yU "° changes by ±0.9 x 10~ 10 from the value for hi = 0. 
Thus, using the LMD±V form factor, our estimate reads 

a L ^° = +5.8 (1.0) x 10~ 10 , (12) 

where the error includes the variation of the parameters and the intrinsic 
model dependence. A similar short-distance analysis in the framework of 
large- Nq QCD and including quark mass corrections for the form factors for 
the j] and 77' was beyond the scope of Re£, We therefore used VMD form 
factors fitted to the relevant CLEO data liil to obtain our final estimate 
a L byL; Ps = fl iw. + a ^ L ;,| vMD + a ™'| VMD = +8.3 (1.2) x 10- 10 . (13) 

We think that an error of 15 % for the pseudoscalar pole contribution is 
reasonable, since we impose many theoretical constraints from long and short 
distances on the form factors. Furthermore, we use experimental information 
whenever available, e.g. the CLEO data for J r 7r o 7 * 7 » (Q 2 , 0) and the decay rate 
7T° — > e + e~. A better measurement of the latter decay could considerably 
reduce the error within the LMD±V ansatz, i.e. the bounds on the parameter 
hi. Our two-dimensional integral representation moreover shows that there 
are no dangerous cancellations, at least in the main contribution aJl byL ' r ^\ 

3 EFT approach to hadronic light-by-light scattering 

In Ref. @ we discussed an EFT approach to <jLbyL ; kad k^g^ on an efFective La- 
grangian_that describes the physics of the Standard Model well below 1 GeV, 
see also E_3. It includes photons, light leptons, and the pseudoscalar mesons and 
obeys chiral symmetn^and U(l) gauge invariance. All pieces of C e g are avail- 
able in the literature Ej and generalizing the chiral power counting by treating 
e, to, and fermion bilinears as order p, one can write £ e g- = +_ + . . . , 
where the relevant terms are given by (for the notation, see Ref. 113) 

£ (2) = - J F^F^ +^{ifl- m)tp + e 2 C{QU+QU) 
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F 2 

+ {{d^U+d^U) +2B (M + U + U + M)) , (14) 
£(4) = F ,u Aa g 0nO + . . . t (15) 

lz7r.ro 

£(6) = 4^b X ^ l5lp + " ' ■ (16) 
The leading contribution to a L 1 byL;had , of order 0(p 6 ), is given by a loop 



of charged pions wit 
value -4.5 x 1CT 10 . 



point-like electromagnetic vertices. It is finite with the 
Since this contribution involves a loop of hadrons, it is 
subleading in the large- Nc expansion E3. 

At leading order in Nc and at order p s in , we encounter the divergent 
pion-pole contribution, diagrams (a) and (b) of Fig. ||, involving two WZW 
vertices from Eq. (fL5"|). Note that the diagram (c) is actually finite. The diver- 
gences of the triangular subgraphs in the diagrams (a) and (b) are removed 
by inserting the counterterm \ from Eq. (|l6|), see the one- loop diagrams (d) 
and (e). Finally, there is an overall divergence of the two- loop diagrams (a) 
and (b) that is removed by a local counterterm denoted by k, diagram (f). 
Note that at order p s there will also be additional contributions, if we re- 
place the point-like vertices in the charged pion-loop by vertices from £' 4 ' , 
e.g. involving the low-energy constants Lg and L w . These contributions will 
be subleading in Nc, although they may be enhanced by logarithms. 




Figure 3. The graphs contributing to a^ byL; " at lowest order in the effective field theory. 

Since the EFT involves a local contribution n to a M , we will not be able 
to give a numerical prediction for a^ byL;had . Nevertheless, we can learn a few 
things by considering the leading logarithms that are in addition enhanced by 
a factor Nc and which can be calculated using the RG tZL The expression of 
the renormalized contribution to arising from the graphs of Fig. || reads 



a 



LbyL;7T 



= H(m/n) + x(v)J(m/n) + K(p). (17) 

The dependence on the subtraction scale fj, in the two-loop function reads 
H(mj[L) = X) P =o l 2 mP ( TO // z )' an d m the one-loop function J(m/fi) = 
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J2 q =o i jq \n- q (m/(i). The physical quantity a^ byUw ° satisfies fj,[da p byU7 ' a /d/j] = 
0. Since the dimcnsionless coefficients h p and j q do not depend on /i, but 
are functions of the ratios M^o/m and F v /m, we have T>h p = and T>j q = 
0, where V = md/dm + M^od/dM^o + F^d/dF^. Finally, the conditions 
T>x(n) — and T>k((i) = hold. We thus obtain the two RG equations 



h 2 = 



dn(fi) 



-7 x io + + hi 



(18) 



with 7 X = n[dx{p) I dp] = Nc, see Ref. 111. The coefficient ji of the logarithm 
in the sum of the one-loop graphs (d) and (e) can easily be worked out 0. In 
this way we obtain from Eq. (En) 



C, C = +3 



in agreement with the result given in Section |[ 



Nc_ 
12tt 



m 



(19) 



4 Further comments and conclusions 



The EFT and large- Nc analysis tells us that we can write 



LbyL;had 




TO 7T± m K ± 
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ln 2 ^ 



, Mo 
ci in — 

m 



Co 
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Mo 



x log's + O 



m 

Tfi ^ 
— j- x Nq x log's 

Mo J 



, (20) 



where u f(m n ± /to, m K ± /m) = —0.038 represents the charged pion and kaon- 
loop that is formally of order one in the chiral and Nc counting. We have 
denoted by fio some hadronic scale, e.g. M p . The coefficient of the log-square 
term in the second line is universal and of order Nc, since F„ = 0(y/Nc)- The 
remaining parts of the coefficient c\ have recently been calculated in Ref. 113 
in the MS-scheme: c\ = — [2x(/io)/3 — 0.904]. The low-energy constant x 
can be determined from the decay 7r° — > e + e~. However, due to the big 
experimental error, this leads to a large uncertainty in x and thus in a^ byLi,u " 1 . 
Assuming lepton universality, one can obtaiiix also from the decay ij — ► /i + /J~, 
with better precision x{M p ) = 1.75^^'qq. E3 One observes that although the 
logarithms are sizeable, \n(M p /m) = 1.98, and give the ballpark of the final 
result, there is some cancellation between the log-square and the log-term in 
a LbyL;had_ ^ ^ e faofc the form factors based on large- Nc QCD that we discussed 
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in Section g give values for the constant Co in Eq. (|2C]) that are of order one, i.e. 
of natural size. Furthermore, in addition to ir° — > e + e~ , we have also taken 
other experimental and theoretical constraints on J r 7r o 7 » 7 » into account. 

It has been argued that a ^ yL ' had has to be positive, because of quark- 
hadron duality, since also the (constituent) quark-loop leads to a positive 
result. We do not agree with this argument for the following reason. Equa- 
tion ( |20| ) tells us that at leading order in Nc any effective theory or model of 
QCD has to show the behavior a ^ uhad ~ (a/ir) 3 Nc[Ncm 2 /(48ir 2 FZ)} In 2 A, 
with a universal coefficient, if pie sends the cutoff A to infinity. From the 
analytical result given in Ref. Ej for the quark-loop one obtains the behav- 
ior a LbyL;c Q M ^ ( a / -k) 3 N c (m 2 / M^) + . . ., for Mq > to, if we interpret the 
constituent quark mass Mq as a hadronic cutoff. Note that the quark-loop is 
also leading in large Nc- lJ Even though one may argue that Nc/16^ 2 F 2 in 
C can be replaced by 1/Mq, the log-square term is not correctly reproduced 
with this model. Therefore, the constituent quark model cannot serve as a 
reliable description for the dominant contribution to a!f yLlhad . On the other 
hand, the coefficient C could, a priori, have any sign, since it is determined 
purely within the EFT, where it is related to the decay amplitude 7r° — > e + e~. 

In conclusion, the pseudoscalar exchange contribution a^yLiPs seems i Q DC 
under control at the 15 % level. In addition, the EFT and large-Ap analysis 
shows a systematic approach to a Jf yL;had and yields the leading and next-to- 
leading logarithmic terms, enhanced by a factor Nc- On the other hand, these 
terms tend to cancel each other to some extent. Furthermore there remains 
the issue of the other contributions to a Lb y L ; had -p or ms tance, the charged pion- 
loop is the leading term in the low-energy expansion, however, after dressing 
the couplings to the photons with form factors, there is a huge suppression, 
by a factor of 3 to 10, although the model-dependent results after the dressing 
differ by about 1 x 10~ 10 . aa Similar statements apply to the dressed quark- 
loop, that appears in these models. Taking all these uncertainties into account 
by adding the model-dependent errors linearly, my conservative estimate for 
the full hadronic light-by-light scattering contribution to a M is as follows: 

a L byL ;had = + g ( 4 ) x 10 -10 _ ( 21 ) 
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